Mean Escape Time in a System with Stochastic Volatility 
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We study the mean escape time in a market model with stochastic volatility. The process followed 
by the volatility is the Cox Ingersoll and Ross process which is widely used to model stock price 
fluctuations. The market model can be considered as a generalization of the Heston model, where 
the geometric Brownian motion is replaced by a random walk in the presence of a cubic nonlinearity. 
We investigate the statistical properties of the escape time of the returns, from a given interval, as 
a function of the three parameters of the model. We find that the noise can have a stabilizing effect 
on the system, as long as the global noise is not too high with respect to the effective potential 
barrier experienced by a fictitious Brownian particle. We compare the probability density function 
of the return escape times of the model with those obtained from real market data. We find that 
they fit very well. 

PACS numbers: 89.65.Gh; 02.50.-r; 05.40.-a; 89.75.-k 



I. INTRODUCTION 

The noise in physical system gives rise to interesting 
and sometimes counterintuitive effects. The stochastic 
resonance and the noise enhanced stability are two ex- 
amples of noise activated phenomena that have been ex- 
tensively studied in a wide variety of natural and phys- 
ical systems such as lasers, spin systems, chemical and 
biological complex systems [HHl Q • Specifically the acti- 
vated escape from a metastable state is important in the 
description of the dynamics of non-equilibrium complex 
systems d, Q • Recently there has been a growing inter- 
est in the application of complex systems methodology 
to model social systems. In particular the application 
of statistical physics for modeling the behavior of finan- 
cial markets has given rise to a new field called econo- 
physics 0,0,1113. 

The stock price evolution is indeed driven by the inter- 
action of a great number of traders. Each one follows his 
own strategy in order to maximize his profit. There are 
fundamental traders who try to invest in solid company, 
speculators who ever try to exploit arbitrage opportunity 
and also noise traders who act in a non-rational way. All 
these considerations allow us to say that the market can 
be thought as a complex system where the rationality 
and the arbitrariness of human decisions are modeled by 
using stochastic processes. 

The price of financial time series was modeled as a ran- 
dom walk, for the first time, by Bachelier His model 
provides only a rough approximation of the real behav- 
ior of the price time series. Indeed it doesn't reproduce 
some of the stylized facts of the financial markets: (i) 
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the distribution of relative price variation (price return) 
has fat tails, showing strong non-Gaussianity [7|,|9|; (ii) 
the standard deviation of the return time series, called 
volatility, is a stochastic process itself characterized by 
long memory and clustering [S [llj ; (iii) autocorrelations 
of asset returns are often negligible [9]. 

A popular model proposed to characterize the stochas- 
tic nature of the volatility is the Heston model [11] , where 
the volatility is described by a process known as the Cox, 
Ingersoll and Ross (CIR) process |13| an d in mathe- 
matical statistics as the Feller process [l4j . The model 
has been recently investigated by econophysicists [H, Gil 
and solved analytically [Tt], [H| . Models of financial mar- 
kets reproducing the most prominent features of statis- 
tical properties of stock market fluctuations and whose 
dynamics is governed by non-linear stochastic differen- 
tial equations have been considered recently in litera- 
ture [lS [H HH [H, HI H, US HI. Moreover finan- 
cial markets present days of normal activity and extreme 
days of crashes and rallies characterized by different be- 
haviors for the volatility. The question whether extreme 
days are outliers or not is still debated. This research 
topic has been addressed both by physicists [26| and 
economists [27l ] . 

A Langevin approach to the market dynamics, where 
market crisis was modeled through the use of a cu- 
bic potential with a metastable state, was already pro- 
posed d, |23|, [2J, [25[ . There feedback effects on the price 
fluctuations were considered in a stochastic dynamical 
equation for instantaneous returns. The evolution inside 
the metastable state represents the normal market behav- 
ior, while the escape from the metastable state represents 
the beginning of a crisis. 

Systems with metastable states are ubiquitous in 
physics. Such systems have been extensively studied. In 
particular it has been proven that the noise can have a 
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stabilizing effect on these systems @, H| . To the best of 
our knowledge all models proposed up to now to study 
the escape from a metastable state contain only a con- 
stant noise intensity, which represents in econophysics 
the volatility. Recently theoretical and empirical inves- 
tigations have been done on the mean exit time (MET) 
of financial time series [H, [2{|, that is the mean time 
when the stochastic process leaves, for the first time, a 
given interval. The authors investigated the MET of as- 
set prices outside a given interval of size L, and they 
found that the MET follows a quadratic growth in terms 
of the region size L. Their theoretical investigation was 
done within the formalism of the continuous time random 
walk. Within the same formalism the statistical proper- 
ties of the waiting times for high-frequency financial data 
have been investigated in refs. [30(. 

In this work we model the volatility with the CIR 
process and investigate the statistical properties of the 
escape times when both an effective potential with a 
metastable state and the CIR stochastic volatility are 
present. Our study provides a natural evolution of the 
models with constant volatility. The analysis has the 
purpose to investigate the role of the noise in financial 
market extending a popular market model, and to pro- 
vide also a starting model for physical systems under the 
influence of a fluctuating noise intensity. The paper is 
organized as follows. In the next section the modified 
Heston model and the noise enhanced stability effect are 
described. In the third section the results for two ex- 
treme cases of this model are reported. In section IV we 
comment the results for the general case and the proba- 
bility density function of the escape time of the returns, 
obtained by our model, is compared with that extracted 
from experimental data of a real market. In the final 
section wc draw our conclusions. 



II. THE MODIFIED HESTON MODEL AND 
THE NES EFFECT 

The Heston model, which describes the dynamics of 
stock prices p(t) as a geometric Brownian motion with 
the volatility given by the CIR mean-reverting process, 
is defined by the following Ito stochastic differential equa- 
tions 

dp(t) = up dt + a(t)p dWx(t) (1) 
dv(t) = a(b-v(t)) dt + c y/v(t?j dW 2 (t), (2) 

where a(t) is the time-dependent volatility, v(t) — a 2 (t) 
and Wi(t) are uncorrelated Wiener processes with the 
usual statistical properties 

< dW t >= 0, < dWi(t)dW 3 (t') >= dt 6{t - t') Si 3 . 

In Eq. (1) fi represents a drift at macroeconomic scales. 
In Eq. ^ the volatility a(t) — y/v(t) reverts towards 



a macroeconomic long time term given by the mean 
squared value b, with a relaxation time a -1 . Here c 
is the amplitude of volatility fluctuations often called 
the volatility of volatility. By introducing log-returns 
x(t) = ln(p(t)/p(0)) in a time window [0, t] and using 
Itd's formula [3l|, we obtain the stochastic differential 
equation (SDE) for x(t) 

dx{t) = Qj, - v(t)/2) dt + y/vjt) dWx{t). (4) 

Here we consider a generalization of the Heston model, 
by replacing the geometric Brownian motion with a ran- 
dom walk in the presence of a cubic nonlinearity. This 
generalization represents a fictitious " Brownian particle" 
moving in an effective potential with a metastable state, 
in order to model those systems with two different dy- 
namical regimes like financial markets in normal activity 
and extreme days (23l . l24l . |25| . The equations of the new 
model are 



dx{t) = - \ — + -^J M + VW) dWx(t) (5) 
dv(t) = a(b-v(t)) dt + cy/vjt) dW 2 (t), (6) 

where U (x) = Ax 3 + Bx 2 is the effective cubic potential 
shown in Fig. [1] with A — 2 and B = 3. From now on we 
call the first and second term of Eq. ^ reverting term 
and noise term respectively. 
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FIG. 1: Cubic potential used in the dynamical equation for 
the variable x(t). The points evidenced in the figure indicate 
the starting positions xo used for the simulations illustrated 
in section UTTl (white circles) and section ITVl (black square). 

Let us call xm the abscissa of the potential maximum 
and Xi the cross point between the potential and the x 
axes. The intervals x < Xi and I — [xi,xm] are clearly 
regions of instability for the system. In systems with 
a metastable state, the random fluctuations can origi- 
nate the noise enhanced stability (NES) phenomenon, an 



3 



interesting effect that increases the stability, enhancing 
the lifetime of the metastable state @, Q . The mean es- 
cape time t for a Brownian particle moving throughout 
a barrier AU is g iven by the the well known exponential 
Kramers law Wl [321 



Cexp 



AU 



(7) 



where r is a monotonically decreasing function of the 
noise intensity v, and C is a prefactor which depends 
on the potential profile. This is true only if the ran- 
dom walk starts from initial positions inside the poten- 
tial well. When the starting position is chosen in the 
instability region x < xm, t exhibits an enhancement 
behavior, with respect to the deterministic escape time, 
as a function of v. This is the NES effect and it can 
be explained by considering the barrier "seen" by the 
Brownian particle starting at the initial position xq, that 
is AUin = U(x max ) — U(xq). Moreover AUi n is less than 
AU as long as the starting position xq lies into the in- 
terval I — [xj,xm]- Therefore for a Brownian particle 
starting from an unstable initial position it is more likely 
to enter into the well than to escape from, once the parti- 
cle has entered. So a small amount of noise can increase 
the lifetime of the metastable state. For a detailed dis- 
cussion on this point and different dynamical regimes see 
Refs. [3|. When the noise intensity v is much greater than 
AU, the Kramers behavior is recovered. 

We have already investigated the statistical properties 
of the escape times for models with stochastic volatility, 
namely the Heston model and the GARCH model [lq . 
We found that the probability density function of the es- 
cape times obtained in the Heston model exhibits a bet- 
ter agreement with real data than that calculated in the 
GARCH model. Here, by considering the modified Hes- 
ton model (Eqs. (5) and ([6|)), characterized by a stochas- 
tic volatility and a nonlinear Langevin equation for the 
returns, we study the mean escape time as a function of 
the model parameters a, b and c. In particular we in- 
vestigate whether it is possible to observe some kind of 
nonmonotonic behavior such that observed for r vs. v 
in the NES effect with constant volatility v. We find a 
nonmonotonic behavior of the mean escape time (MET) , 
r, as a function of the model parameters. Within this 
behaviour we recognize the enhancement of t as a NES 
effect in a broad sense. 

Our modified Heston model has two limit regimes, cor- 
responding to the cases a — 0, with only the noise term 
in the equation for the volatility v(t), and c = with only 
the reverting term in the same equation. This last case 
corresponds to the usual parametric constant volatility 
regime. In fact, apart from an exponential transient, the 
volatility reaches the asymptotic value b. The NES effect 
should be observable in the latter case as a function of 
b, which is the average volatility. In fact, in this case we 
recover the motion of a Brownian particle in a fixed cubic 
potential with a metastable state and an enhancement of 



its lifetime for particular initial conditions. 

For this purpose we perform simulations by integrat- 
ing numerically the equations (5) and ([6]) using a time 
step At = 0.01, and, as potential parameters, A = 2.0 
and B = 3.0. All the integration steps yielding negative 
values for v were rejected and repeated. The simulations 
were performed placing the walker in the initial positions 
xq located in the unstable region [xi, xm] (see Fig.[l} and 
using an absorbing barrier at x — —6.0. Each time the 
walker hits the barrier, the escape time is registered and 
another simulation starts, placing the walker at the same 
initial position xq, but using the volatility value of the 
barrier hitting time. 



III. LIMIT CASES 

First of all we present the result obtained in the limit 
cases where we have only one of the two terms in the 
CIR equation ([6]). Namely: (a) only the reverting term 
(o/0,i/0,c=0, revert-only case), and (b) only the 
noise term (a — 0, b whatever, c ^ 0, noise-only case) are 
present. In the case (a) the volatility is practically con- 
stant and equal to b, apart from an exponential transient 
that is negligible for times t ^S> a -1 . 

The mean escape time as a function of b is plotted in 
Fig. for the seven different initial positions indicated 
in Fig. [T] (white circles in that figure). The curves are 
averaged over 10 5 escape events. The nonmonotonic be- 
havior is present. The escape time increases by increas- 
ing b until it reaches a maximum. After the maximum, 
when the values of b are much greater than the potential 
barrier height, the Kramers behavior is recovered. The 
nonmonotonic behavior is more evident for starting po- 
sitions near the maximum of the potential. For starting 
positions xo lying in the interval I — [xj,Xm], the ini- 
tial plateau is an artifact of the calculus. Indeed in the 
theoretical solution for the constant volatility case, r di- 
verges as the noise intensity approaches zero [H, 0] • With 
such a noise intensity the escape from the well is a very 
unlikely event. We should require an infinite number of 
simulation steps in order to observe an event pushing the 
particle into the well. 

In the case (b) , v is a multiplicative stochastic process 
with standard deviation equal to c. The random term in 
Eq. ^ is given by the product of the white noise and the 
square root of the process itself v. In Fig. [2}d we plot r as 
a function of the parameter c. The nonmonotonic behav- 
ior is absent and the MET decreases monotonically, but 
differently from the Kramers behavior. The curves are 
averaged over 10 s escape events. In this case a greater 
number of events is required to eliminate the fluctuations 
present in the curves. Once again we cannot draw any 
conclusions for the constant behavior of t, when the val- 
ues of the parameter c are very small, because of the finite 
number of simulation steps [3j . It is worthwhile pointing 
out that the NES effect is not observable as a function 
of the parameter c, if a = 0. Therefore, the presence of 
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FIG. 2: Mean escape time r for the seven different start- 
ing positions of Fig. [T] (white circles in that figure). Limit 
cases where only one of the two terms in the CIR equation is 
present: (a) r vs. b, when only the reverting term is present 
(a — 10~ 2 , c = 0), and (b) r vs. c, when only the noise term 
is present (a = 0, b = 10 -2 ). The different starting posi- 
tions xo from top to bottom are: —1.1, —1.2, —1.3, —1.4, 
— 1.5, —1.55, —1.60. The parameters b and c are dimension- 
less, while a, measured in arbitrary units, has dimension of 
the inverse time. Inset: detail of the curve with xo = —1.10. 



the reverting term affects the behavior of r in the do- 
main of the noise term of the volatility and it regulates 
the transition from nonmonotonic to monotonic regimes 
of MET. Moreover in this noise-only regime the volatility 
is proportional to the square of the Wiener process and 
therefore the fluctuations during the monotonic decreas- 
ing behavior of r are very large. We see three different 
dynamical regimes: (i) for low values of the parameter 
c, with respect to the height of the potential barrier, we 
have a constant behavior, (ii) for intermediate values of 
c we have fluctuations with decreasing monotonic behav- 
ior, (iii) for values of c greater than 1, which is exactly 



the height of the barrier, we get very small and constant 
values of t close to zero. This behavior is mainly due to 
the presence of the Ito term in Eq. (5) for the log-returns 
x(t). In fact, because of this term we have a fluctuat- 
ing potential, obtained by the previous one by adding 
a linear term: U — U(x) + (xv)/2. The effect of the 
positive linear term is to modify randomly the potential 
shape in such a way that the potential barrier disappears 
for greater values of the volatility v(t). Specifically for 
values of v < 3 the potential barrier is always present, 
while for v > 3 it disappears. This produces a random 
enhancement of the escape process with a consequent de- 
creasing behavior of the average escape time. 



IV. GENERIC CASE 

In order to present our results for the generic case, 
where both the reverting and noise terms of the CIR 
equation are present, we focus on a single starting posi- 
tion, xq = —1.25, which is located in the middle of the 
interval [xi,xm]- We analyze the escape time through 
the barrier using different values for parameters a, b and 
c. We note that the average escape time r is measured 
in arbitrary units (a. u.), the parameter a (measured in 
a.u. too) has the dimension of the inverse time, while b 
and c are dimensionless. 

As a first result we present the behavior observed for 
t as a function of the reverting level b. In Fig. [3] we 
show the curves averaged over 10 5 escape events. Each 
panel corresponds to a different value of c. Inside each 
panel different curves correspond to different values of a 
spanning seven orders of magnitude. The nonmonotonic 
shape, characteristic of the NES effect, is clearly shown in 
Fig. This behavior is shifted towards higher values of 
b as the parameter a decreases, and it is always present. 
In Fig. [3J;, which corresponds to a much greater value of 
c (c = 30), all the curves are monotonic but with a large 
plateau. So an increase in the value of c causes the NES 
effect to disappear. 

To understand this behavior let us note that the pa- 
rameters a and c play a regulatory role in Eq. For 
a c the drift term is predominant while for a <C c 
the dynamics is driven by the noise term, unless the pa- 
rameter b takes great values. In fact in Fig. the non- 
monotonic behavior is observed for a -C c, provided that 
b^> c. For increasing values of a the system approaches 
the revert-only regime and we recover the behavior shown 
in Fig. [2^. For a <C c the shape of the curves changes. 
The mean escape time r is almost constant, and only for 
very high values of b we observe a decreasing behavior. 
This happens because for smaller values of a the revert- 
ing term becomes negligible in comparison with the noise 
term and the dependence on b becomes weaker. Only 
when b is large enough the reverting term assumes val- 
ues that are no more negligible with respect to the noise 
term and we can observe again a dependence of r on b. 
By increasing the value of c we observe the nonmonotonic 
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FIG. 3: Mean escape time r as a function of reverting level b. 
Each panel corresponds to a different value of c, specifically 
(a) c = 10~ 2 , (b) c = 10 and (c) c = 30. Inside each panel 
different curves correspond to the following values of a: (a) 
black circle 10 -7 , white circle 10~ 6 , black triangle down 10~ 5 , 
white triangle up 10 -4 , black square 10 -3 , white square 10 -2 , 
black diamond 10 _1 ; (b) and (c) black circle 10 -4 , white circle 
10 _1 , black triangle down 1, white triangle up 10, black square 
10 2 . 



behavior only for a very great value of the parameter a, 
that is for a = 100 » c (see Fig. [3b). For further increase 
of parameter c (see Fig. [5b), the noise experienced by the 
system is much greater than the effective potential bar- 
rier "seen" by the fictitious Brownian particle and the 
NES effect is never observable. Moreover we are near a 
noise-only regime, and we can say that the magnitude of 
c is so high as to saturate the system. 

In summary the NES effect can be observed as a func- 
tion of the volatility reverting level b, the effect being 
modulated by the parameter (ab) / c. The phenomenon 
disappears if the noise term is predominant in compar- 
ison with the reverting term. Moreover the effect is no 
more observable if the parameter c pushes the system 
towards a too noisy region. 

As a second step we study the dependence of r on 
the noise intensity c. Fig. 3] shows the curves of r vs c, 
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FIG. 4: Mean escape time r as a function of the noise intensity 
c. Each panel corresponds to a different value of b, specifically 
(a) b — 10 -2 and (b) b = 10. Inside each panel different 
curves correspond to the following values of a: black circle 
10" 3 , white circle 3 x 10~ 3 , black triangle down 9 x 10" 3 , 
white triangle up 2.7 x 10~ 2 , black square 8.1 x 10 -2 , white 
square 2.4 x 10 _1 , black diamond 7.3 x 10 _1 , white diamond 
2.2, black triangle up 6.6. 
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averaged over 10 5 escape events. Each panel corresponds 
to a different value of b. Inside each panel different curves 
correspond to different values of a. The shape of the 
curves is similar to that observed in Fig. [31 For small b 
(panel a) we observe a nonmonotonic behavior, while for 
great b (panel b) the curves are monotonic but with a 
large plateau. Let us recall the results of Fig. there 
was no NES effect in the noise-only case. So for small 
values of a, when the reverting term is negligible, the 
absence of the nonmonotonic behavior is expected. By 
increasing a the nonmonotonic behavior is recovered (see 
Fig. Once again, if one of the parameter pushes 

the system into a high noise region, the nonmonotonic 
behavior disappears (see Fig. [3b). Specifically if b is high, 
the reverting term drives the system towards values of 
volatility that are outside the region where the NES effect 
is observable. Indeed a direct inspection of Fig. [3] shows 
that the value of b used in Fig. 2Jd is located after the 
maximum of r for all values of a and c. 

In summary when the noise term is coupled to the 
reverting term we observe the NES effect on the variable 
c. The effect disappears if b is so high as to saturate the 
system. As a last result we discuss the behavior observed 
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FIG. 5: Mean escape time r as a function of the reverting 
rate a for different values of b and c. Curves with the same 
color correspond to the same value of c, while curves with the 
same symbol correspond to the same value of 6. Specifically 
the values of b used are: 10 -3 circle, 1 triangle, 10 2 square . 
The values of c used are: 10" 2 black, 10" 1 white, 1 gray. 



for r as a function of the reverting rate a. This allow us 
to observe the transition between the two regimes of the 
process discussed above: the noise-only regime and the 
revert-only regime. The results are reported in Fig. [5] for 
three different values of b and three different values of 
c. To reduce the fluctuations in all the curves when the 
parameter a becomes small, we performed simulations 
by averaging on 10 6 escape events. It is worthwhile to 
note that for values of the parameter a < 10~ 5 , we enter 
in the noise-only regime, which characterizes one of the 
limit cases discussed in section III. Curves with the same 



color correspond to the same value of c while curves with 
the same symbol correspond to the same value of b. 

The system tends to the noise-only regime for lower 
values of a and to the revert-only regime for higher values 
of a. On the right end of Fig. \E\ the curves correspond- 
ing to the same value of b tend to group together. The 
values of r the curves approach reflect the nonmonotonic 
behavior observed in Fig. [3Ji. Indeed all the curves corre- 
sponding to the intermediate value of & (b = 1) approach 
a value of r, which almost corresponds to the maximum 
value of MET in Fig. [3^ (we note that the behavior for 
a = 10 _1 coincides with that for a = 1, even if it is not 
reported in Fig. [3^) . This value of r is greater than that 
reached by the curves corresponding to other two, lower 
and greater, values of & (10 -3 and 10 2 respectively). Con- 
versely on the left end the curves corresponding to the 
same value of c tend to group together. It is worth noting 
that in this last case the curves with the highest value 
of c, namely c = 1, show greater fluctuations as those 
observed in all the previous cases where the noise term 
is predominant (see for example Fig. [2b)- 

It is interesting to show, for our model (Eqs. (5) 
and (6)), some of the well-established statistical signa- 
tures of the financial time series, such as the probability 
density function (PDF) of the stock price returns, the 
PDF of the volatility and the return correlation. In Fig.[B] 
we show the PDF of the returns. To characterize quan- 
titatively this PDF with regard to the width, the asym- 
metry and the fatness of the distribution, we calculate 
the mean value < Ax >, the variance (JAxj the skewness 
K3, and the kurtosis We obtain the following val- 
ues: < Ax >= -0.162, cr Ax = 0.348, n 3 = -1.958, 
K4 = 5.374. These statistical quantities clearly show 
the asymmetry of the distribution and its leptokurtic na- 
ture observed in real market data. In fact, the empirical 
PDF is characterized by a narrow and large maximum, 
and fat tails in comparison with the Gaussian distribu- 
tion 0,1- Specifically we note that the value of the kur- 
tosis K4 = 5.374, which gives a measure of the distance 
between our distribution and a Gaussian one, is of the 
same order of magnitude of that obtained for S&cP 500 
for daily prices (see Fig. 7. 2 on page 114 of Ref. '8]). The 
presence of the asymmetry is very interesting and it will 
be subject of future investigations. This asymmetry is 
due to the nonlinearity introduced in the model through 
the cubic potential (see Fig. [T|) . Of course a comparison 
between the PDF of real data and that obtained from 
the model requires further investigations on the dynam- 
ical behavior of the system, as a function of the model 
parameters. In the following Fig. [7] we show the PDF 
of the volatility for our model, and we can see a log- 
normal behavior as that observed approximately in real 
market data. Finally in Fig. [5] we show the correlation 
function of the returns. As we can see the autocorrela- 
tions of the asset returns are insignificant, except for very 
small time scale for which microstructure effects come 
into play. This is in agreement with one of the stylized 
empirical facts emerging from the statistical analysis of 
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FIG. 6: Probability density function of the stock price returns 
for our model (Eqs. (5) and (6)). The values of the param- 
eters are: a = 10 _1 , b = 4.5, c = 2 x 10 _1 . The potential 
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FIG. 8: Correlation function of the returns for our model 
(Eqs. (5) and (6)). The values of the parameters are the 
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FIG. 7: Probability density function of the volatility for our 
model (Eqs. (5) and (6)). The values of the parameters are 
the same of Fig. [6] Inset: semilog plot of the PDF of volatility 
in a longer time scale. 
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FIG. 9: Probability density function of the escape time of 
the returns from simulation (solid line), and from real data 
(black circle). The values of the parameters are: Axi — 
— O.Ictax, Axj = — 1-Octax, x = -1.25, x abs = -6.0, 
vstart = 8.62 x 1(T 5 , a = 10" 1 , b = 4.5, c = 2 x HT 1 . 
The potential parameters are: A — 2 and B = 3. 



price variations in various types of financial markets [9(. 
A quantitative agreement of the PDF of volatility and 
the correlation of returns with the corresponding quanti- 
ties obtained from real market data is subject of further 
studies. 

Our last investigation concerns the PDF of the escape 
time of the returns, which is the main focus of our paper. 
By using our model (Eqs. (5) and ©), we calculate the 
probability density function for the escape time of the 
returns. We define two thresholds, Axi and Ace/, which 
represent the start point and the end point for calculat- 
ing r respectively. When the return series reaches the 
value Axi, the simulation starts to count the time r and 
it stops when the threshold Ax j is crossed. In order to fix 
the values of the two thresholds we consider the standard 
deviation (SD) a& x of the return series over a long time 



period corresponding to that of the real data. Specifically 
a Ax is the average of the standard deviations a n observed 
for each stock during the above mentioned whole time pe- 
riod (n is the stock index, varying between 1 and 1071). 
Then we set Axi = —0.1 <j& x and Axf = —1.0 <j& x - 
We perform our simulations obtaining a number of time 
series of the returns equal to the number of stocks con- 
sidered, which is 1071. The initial position is xq = —1.25 
and the absorbing barrier is at x a bs — —6.0. For the 
CIR stochastic process v, we choose v sta rt — 8.62 x 10~ 5 , 
a = lQ-\ b = 4.5 and c = 2 x 10 -1 . The choice of this 
parameter data set is not based on a fitting procedure 
as that used for example in Ref. [l|j]. There the mini- 
mization of the mean square deviation between the PDF 
of the returns, extracted from financial data, and that 
obtained theoretically is done. We choose the parameter 
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set in the range in which we observe the nonmonotonic 
behaviour of the mean escape time of the price returns as 
a function of the parameters b and c. Then by a trial and 
error procedure we select the values of the parameters a, 
b, and c for which we obtain the best fitting between the 
PDF of the escape times calculated from the modified 
Heston model (Eqs. (5) and (6)) and that obtained from 
the time series of real market data. We report the re- 
sults in Fig. [9] Of course a better quantitative fitting- 
procedure could be done, by considering also the poten- 
tial parameters. This detailed analysis will be done in a 
forthcoming paper. 

As real data we use the daily closure prices for 1071 
stocks traded at the NYSE and continuously present in 
the 12-year period 1987 - 1998 (3030 trading days). 
The same data set was used in previous investigations 
Il7l . l3-'j| , |34| . From this data set we obtain the time 
series of the returns and we calculate the time to hit a 
fixed threshold starting from a fixed initial position. The 
two thresholds were chosen as a fraction of the average 
standard deviation a^ x on the whole time period, as we 
have done in simulations. The agreement between real 
data and those obtained from our model is very good. 
We note that at high escape times the statistical accu- 
racy is worse because of few data with high values. The 
parameter values of the CIR process for which we ob- 
tain good agreement between real and theoretical data 
are in the range in which we observe the nonmonotonic 
behavior of MET (see Fig. [3ja,). This means that in this 
parameter region we observe a stabilizing effect of the 
noise on the prices in the time windows for which we 
have a variation of returns between the two fixed values 
Axi and Ax f . This encourages us to extend our analysis 
to large amounts of financial data and to explore other 
parameter regions of the model. 



V. CONCLUSIONS 

We studied the mean escape time in a market model 
with a cubic nonlincarity coupled with a stochastic 
volatility described by the Cox-Ingersoll-Ross equation. 
In the CIR process the volatility has fluctuations of in- 
tensity c and it reverts to a mean level b at rate a. 

Our results show that as long as the mean level a is 
different from zero it is possible to observe a nonmono- 
tonic behavior of MET as a function of the two model 
parameters b and c. The parameter a regulates the tran- 
sition from a noise-only regime, where reverting term is 
absent or negligible, to a revert-only regime, where the 
noise term is absent or negligible. In the former case, 
the enhancement of MET with a nonmonotonic behavior 
as a function of the model parameters, that is the NES 
effect, is not observable. The curves have a monotonic 
shape with a plateau. Moreover, if one of the parame- 
ters is so big to push the system into a region where the 
noise is greater than the barrier height of the effective 
potential, the effect is no more observable at all. In the 
revert-only regime, instead, the NES phenomenon is re- 
covered. With its regulatory effect, the reverting rate a 
can be used to modulate the intensity of the stabilizing 
effect of the noise observed by varying b and c. In this 
parameter region the probability density function of the 
escape times of the returns fits very well that obtained 
from the experimental data extracted by real market. 
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